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Abstract 

We consider renormalizable SO (10) Yukawa interactions and put the three fermi- 
onic 16-plets into the 3-dimensional irreducible A4 representation. Scanning the 
possible A4 representation assignments to the scalars, we find a unique case which al- 
lows to accommodate the down-quark and charged-lepton masses. Assuming type II 
seesaw dominance, we obtain a viable scenario with the Zee-Wolfenstein neutrino 
mass matrix, i.e., the Majorana mass matrix with a vanishing diagonal. Contribu- 
tions from the charged-lepton mass matrix resolve the well-known problems with 
lepton mixing arising from the vanishing diagonal. In our scenario, fermion masses 
and mixings are well reproduced for both normal and inverted neutrino mass spec- 
tra, and b-T Yukawa unification and definite predictions for the effective mass in 
neutrinoless double-/? decay are obtained. 
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1 Introduction 



Grand Unified Theories (GUTs) based on the gauge group SO (10) pQ are a framework 
for attempts to understand the observed fermion masses and mixings. These theories fea- 
ture a 16-dimensional irreducible representation (irrep), the spinor representation, which 
naturally accommodates all chiral fermions of one Standard Model (SM) generation plus 
a right-handed neutrino. Furthermore, 5*0(10) GUTs allow for type I [2] and type II [3] 
seesaw mechanisms (see also [I]) for explaining the smallness of the light neutrino masses. 

By employing only one scalar in the 10 and one in the 126 irrep of 5*0(10) in 
renormalizable Yukawa couplings, the so-called "Minimal Supersymmetric 5O(10) GUT" 
(MSGUT) [5] has been successful in accounting for all fermion masses and mixings, if 
one focuses solely on its fermion mass matrices. Moreover, this model has built-in the 
gauge-coupling unification of the Minimal Supersymmetric Standard Model (MSSM). In- 
depth studies of the MSGUT have been performed [HI Ej, also with inclusion of small [5] 
and prominent [H] effects of the 120 irrep. Despite its success in reproducing the known 
fermion masses and mixings, it should be stressed, however, that the MSGUT considered 
as a whole is too constrained, as its scalar sector does not admit the vacuum expectation 
values (VEVs) required for the fit in the fermionic sector — see [10] and [7j. 

However, one is not really satisfied by just reproducing known fermion masses and 
mixings, one would also like to explain, for instance, the threefold replication of fermion 
generations, or the peculiar mixing properties of the lepton sector [TT], that is, maxi- 
mal atmospheric and large but non- maximal solar mixing and a small mixing angle #13. 
5O(10) models have not yet been successful in explaining such features, but one may try 
to meet such challenges by considering the possibility of an underlying flavor symmetry 
group G. In GUT models based on the product group 5O(10) x G, the three known 
fermion generations can be assigned to representations of the flavor group G. Our choice 
of G is guided by the following considerations. Since all irreps of abelian groups are one- 
dimensional, only non-abelian groups are suited to explain the existence of more than 
one fermion generation. Furthermore, unlike continuous symmetries, the break-down of 
discrete global symmetries does, in general, not give rise to undesired Goldstone bosons. 
This suggests to stick to non-abelian, discrete flavor groups. Only a few 5O(10) x G mod- 
els, with G being non-abelian and discrete, have been studied so far. For instance, models 
with 5O(10) x 54 [12] and 5O(10) x A4 [13] symmetries have already been investigated. 
In particular, models employing an A4 [T^l [15] flavor symmetry may give tri-bimaximal 
leptonic mixing [16j, but, in general, in such models right- and left-handed fermion fields 
transform differently under A4. However, in 5O(10) GUTs right- and left-handed fermion 
fields have to transform in the same way under [17], since there all chiral matter 
fields of one generation belong to the same 5O(10) irrep. In [13] a non-supersymmetric 
5O(10) x A4 model with type I seesaw dominance has been analyzed, which successfully 
preserves tri-bimaximal leptonic mixing and can accommodate all known fermion masses. 
The quark mixing angles, however, are assumed to be zero. 

In this paper, we investigate the fermionic sector of renormalizable 5O(10) x A4 GUT 
scenarios, with the three fermion families in the three-dimensional irrep of A4, while for 
the 5O(10) scalar irreps occurring in Yukawa couplings we allow all possible A4 irreps. 
We do not discuss the difficult problem of vacuum alignment, but rather assume that we 
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can dispose of the VEVs according to our needs. With this assumption, we will see that 
the SO(10) x A 4 structure enforces the building blocks of the fermion mass matrices to 
consist of diagonal and off-diagonal matrices. A crucial role will play the mass matrices 
of the down quarks and the charged leptons. Requiring solely that the scenario is able 
to reproduce down-quark masses and charged-lepton masses, singles out a unique case 
with respect to the transformation of the scalars under SO (10) x A*. In that unique 
viable scenario, under the assumption of type II seesaw dominance, we will find the 
Zee-Wolfenstein form [181 EH] of the mass matrix of light neutrinos. The well-known 
phenomenological problems [20] of this mass matrix turn out to be completely resolvable 
by contributions to lepton mixing from the charged-lepton sector. We want to stress, 
however, that our usage of the flavor symmetry does not enforce tri-bimaximal mixing 
in the lepton sector. 

Though we have in mind a supersymmetrized scenario, supersymmetry enters our con- 
siderations only via the fermion masses. For the numerics we use masses at the GUT scale, 
which have been obtained through the renormalization group equations of the MSSM. 

Our paper is organized as follows. In Section [2] we summarize the properties of the 
Zee-Wolfenstein neutrino mass matrix. The SO (10) x A4 GUT scenario is developed in 
Section [3j The methods and results of our numerical analysis are discussed in Section @J 
Section [5] is devoted to the conclusions. 

2 The Zee— Wolfenstein mass matrix in a nutshell 

The Zee model generates Majorana neutrino masses at the one- loop level [TBI [21] . Its 
neutrino mass matrix has, in general, non-zero elements on the diagonal. However, with 
a suitable symmetry one can enforce a vanishing diagonal in M. v at the one- loop 
level [19J. This Zee-Wolfenstein neutrino mass matrix is a symmetric 3x3 matrix whose 
diagonal elements are zero: 



From now on we will discuss only this case. The restricted Zee model has the property 
that one can make a basis transformation such that the charged-lepton mass matrix is 
diagonal but the form of M. v given by (pQ) persists. Thus, without loss of generality, we 
will assume in this section that the charged-lepton mass matrix is diagonal. 

In diagonalizing the matrix (fTl) , one can first remove the phases of M. u . These phases 
can be absorbed into the charged-lepton fields. Thus we take the matrix entries a, b, c to 
be real. Since the Zee-Wolfenstein mass matrix is traceless and symmetric one has [22] 



where the Aj denote the real eigenvalues of Ai v . Writing m ; for the masses of the light 
neutrinos, we have rrii = |Aj|. 

In the case of inverted ordering 777,3 < mi < 1712 of the neutrino masses (Am|, = 
ml— ml, Am^ tm = rn| — it has been pointed out in Ref. [20] that the mass matrix ([!]) 
together with Am% <C Am^ tm leads, for all practical purposes, to maximal solar mixing 




(1) 



Ai + A 2 + A 3 = 0, 



(2) 
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#12 = vr/4 and #13 = 0. Furthermore, the atmospheric mixing angle #23 can be chosen to 
be maximal. While the latter two properties are most welcome, maximal solar mixing is 
excluded by more than 5a by experimental data [23J. In [21] it is has been shown that 
deviations from maximal solar mixing are severely constrained through 

, / 1 Am?, , v 

cos20 12 < 7-r— 3 
4 Ami 



'atm 



The neutrino masses are approximately given by 

1 Am 2 



n ' :i - « Ti "n y m 1 ~m 2 ~ y/ A"»L . (4) 
V Am atm 



Thus one obtains 7713 <C mi ~ m 2 and the resulting neutrino mass spectrum exhibits an 
inverted hierarchy. 

In the case of the normal ordering mi < m 2 < m 3 of the neutrino masses (Am?, = 
m\ — m 2 , Am 2 tm = m\ — mf) things are even worse. Although it is now possible to have 
maximal atmospheric mixing and, at the same time, allowing the solar mixing angle to be 
in perfect agreement with experimental data, the mixing angle # 13 turns out to be much 
too large [20] : 

sin 2 13 ~ - . (5) 
The neutrino mass spectrum can be estimated by 



1 /AmL . . 

mi~m 2 ~-m 3 ~Y — - — . (6) 

Therefore, all neutrino masses will be of the same order of magnitude, but the mass 
spectrum cannot be quasi-degenerate. 

Concerning neutrinoless double-/? decay, the relevant observable is the effective Majo- 
rana neutrino mass \(m^p)\ = E»^ei m *l> where U denotes the unitary leptonic mixing 
(PMNS) matrix. The mass \{mpp)\ is equal to the modulus of the (e, e) matrix element 
of Ai u , which is exactly zero in the Zee-Wolfenstein case. Thus the model prevents 
neutrinoless double-/? decay. 

In summary, the Zee-Wolfenstein model is not viable because it does not give a con- 
sistent explanation of all current experimental data of the neutrino sector (i.e. two mass- 
squared differences plus three mixings angles). It is the purpose of this paper to embed 
the Zee-Wolfenstein neutrino mass matrix in an SO (10) GUT. In such an environment, 
the zeros in the diagonal of M. v are not stable under a basis change such that the charged- 
lepton mass matrix becomes diagonal. Therefore, as we will show, contributions from the 
charged-lepton sector can provide the necessary remedy for correcting the too large mixing 
angle #12 in the case of inverted hierarchy and #13 in the case of normal hierarchy [22] . As 
an additional bonus, a non-vanishing | (mpp) \ and, therefore, neutrinoless double-/? decay 
becomes possible. 
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3 The 50(10) x A 4 model 

The tensor product of the SO (10) spinor representation of the fermions is given by [261I2T] 

16® 16 = (10©126) S ©120 AS , (7) 

where the subscripts S and AS refer to symmetric and antisymmetric Yukawa coupling 
matrices, respectively. Renormalizable SO (10) GUTs can generate fermion masses at the 
tree level only by the scalar irreps 10, 120 and 126. 

The 12-element group A4 is popular as a family symmetry in model building — see [H] 
for a selection of the vast A 4 literature and [15] for a review on the group A4 and models. 
It has three one-dimensional irreps and one three-dimensional irrep. The tensor product 
3 © 3 contains all one-dimensional irreps exactly once, but the three-dimensional irrep 
is contained twice. While the Yukawa couplings corresponding to the one-dimensional 
irreps are diagonal and, therefore, symmetric, the couplings of the 3 © 3 G 3 © 3 are 
off-diagonal, but no special symmetry property is fixed. However, Eq. ([7j) suggests to 
choose one three-dimensional irrep with symmetric and the other one with antisymmetric 
tensor indices: 

3 © 3 = (1 © 1' © 1" © 3) s © 3 AS - (8) 

Now we consider 50(10) x A4 and investigate possible Yukawa couplings and fermion 
mass matrices under the assumption that the fermions transform as 16 © 3, which is 
is clearly the only reasonable choice if we want to take advantage of the non-abelian 
character of A4. Equations ([7]) and (|HJ) dictate that the 120 can only transform as a 3 
under A4, while for the for 10 and 126 singlet and triplet irreps of A^ are possible. Let 
us consider the case where the scalars responsible for Yukawa couplings transform as 

10©(1©1' ©1") and 126 ©3. (9) 

Then, in a symbolic way, writing down only the A^ part, the Yukawa couplings are given 
by 

3 3 

h J2 ^ 1)(a ' 1] 16 a 16 a 10, + (10) 

i=l a=l 

(16 2 16 3 + 16 3 16 2 ) T26! + (16 3 16! + 16x16a) 126 2 + (16!l6 2 + 16 2 16i) 126 3 , (11) 

where a is a family index and to = (— 1 + iy/3)/2. Furthermore, we make two assumptions: 

i) All VEVs which occur in the scalars can have independent values. 

ii) Type II seesaw dominates in the neutrino mass matrix. 

These assumptions together with Eq. (Q define the scenario we will investigate in the 
following. 

We furthermore assume that our models can be extended in a suitable way to solve the 
doublet-triplet splitting problem^ Moreover, since we have in mind the MSSM, with only 

1 For instance, the Dimopoulous-Wilczek mechanism [29] and the missing partner mechanism |30j 
provide viable solutions of the doublet-triplet splitting problem in SO(10) GUTs. 
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two Higgs doublets, as the low-energy limit of our SO (10) models we must assume a suit- 
able doublet-doublet splitting as well, which is usually achieved by finetuning [31] . These 
assumptions are not innate to the models presented here but are well-known problems in 
GUTs. 

Let us now derive some consequences of our scenario. Because of assumption i), the 
Yukawa couplings ( flOi) produce diagonal mass terms with three independent entries, one 
for the up-quark mass matrix (q = u) and another one for the down quark mass matrix 
(q = d): 

diag (hiv ql + h 2 v q2 + h 3 v q3 , h x v ql + uh 2 v q2 + uj 2 h 3 v q3 , h t v ql + u?h 2 v q2 + ujh 3 v q3 ) , (12) 

where the v q i are the VEVs appearing in the scalar 10-plets. Next we consider the Yukawa 
couplings ( fTTT) . Again because of assumption i), this Yukawa interaction generates two 
independent off- diagonal contributions to the mass matrices of up and down quarks. 

Studying the system of mass matrices, we find that assumptions i) and ii) lead to 
a decoupling of the up-quark mass matrix M u from the rest of the system. This is so 
because the quark mass matrices are given by M u = H' + F' and Md = H + F, where H' 
and H are independent diagonal matrices, while F' and F are independent off-diagonal 
matrices. Therefore, M u would only be related to the system of mass matrices through 
the neutrino Dirac-mass matrix Md = H' — 3 F', but this relationship is irrelevant due 
to assumption ii). Since M u is a general symmetric matrix independent of the rest of the 
system of mass matrices, the CKM matrix can always be reproduced. The other side of 
the coin is that our scenario loses predictivity because it is neither restricted by the values 
of the up-quark masses nor by the experimental information on the CKM matrix. 

The remaining system of mass matrices which we want to study consists of the mass 
matrices of down-type quarks and charged-leptons, given by 

M d = H + F and M e = H - 3 F, (13) 

respectively, where H is diagonal, while F is off-diagonal, and of the neutrino mass matrix 
M. v of Eq. ([I]). Without loss of generality, H can be assumed to be real, but F and 
A4 U have complex entries^] Note that in view of assumption i) the entries in A4 U are 
independent of F, but and Ai u are coupled via the PMNS matrix 

U = lfjU u with UjMfUt = diag (m e , m M , m T ), M. U U U = diag (mi, m 2 , m 3 ). (14) 

Counting the number of parameters, we find nine absolute values and five phasesjf) while 
the number of observables to be fitted is 11: three charged-lepton masses, three down- 
quark masses, two neutrino mass-squared differences and three lepton mixing angles. The 
fitting procedure and predictions of our scenario will be exposed in the next section. 

We note that the family symmetry has the effect of generating independent diagonal 
and off-diagonal contributions to the quark and lepton mass matrices. Adhering to the 
two assumptions presented above but using other A 4 representations than those of Eq. Q, 
we can find several other scenarios. E.g., with 120-plets, antisymmetric off-diagonal mass 

2 Now we cannot absorb the phases of M. v into the charged-lepton fields since Me is not diagonal. 
3 Of the three phases in M v one can be removed. 
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Leptons 




n +0.0003 

_o,0003 




75.6715+0-^8 

— U.UoUl 


m T 


1292.2 




(7.9 ±0.3) x 10~ 5 


Am ltm 


(2.50 t° f 5 ) x 10- 3 


s 2 


0.31 ±0.025 


S 2 
6 23 


0.50 ±0.065 


S 2 


< 0.0155 



Quarks 


m d 


1.03 ±0.41 




19.6 ±5.2 


m b 





Table 1: Input data (central values and la errors) at the GUT scale of Mqut = 2 x 
10 16 GeV for tan/3 = 10. The charged-fermion masses are taken from [32], except for the 
values of and m s ; these were obtained by taking their low-energy values from [33] and 
scaling them to M GU t- As for Am 2 tm , we use the value obtained in [23] • We have copied 
the remaining input from Table I in [7J. Charged-fermion masses are in units of MeV, 
neutrino mass-squared differences in eV 2 . We have used the abbreviations s\ 2 = sin 2 #12, 
etc. The angles in the left table refer to the PMNS matrix. 

matrix contributions are generated. A list of such cases is presented in Table El There 
we confine ourselves to a maximum of three scalars per 50(10) irrep, the 126-plet must 
always be present to allow a viable type II seesaw neutrino mass matrix and the 10 and 
120 are not present at the same time; the latter condition is for avoiding a proliferation 
of parameters. However, it will turn out that the only viable scenario is the one defined 
via Eq. ©. 

4 The numerical analysis 

We perform a global \ 2 analysis of the S'O(IO) x A4 scenario defined by Eq. (jUJ) and 
assumptions i) and ii) by employing the downhill simplex method [28] . In Table [1] the 
observable quantities Oi are specified in the form 

Oi = Oi ± <r i} (15) 

where Oi and cij denote central values and la deviations, respectively. The index % = 
1, . . . , 11 labels the different observables given in Table [lj The masses in that table refer 
to the mass values at a GUT scale of 2 x 10 16 GeV, obtained via the renormalization 
group equations of the MSSM, for the ratio of Higgs doublet VEVs tan (3 = 100 Writing 

Since we do not have quasi-degenerate neutrino mass spectra, the effect of the renormalization group 
running on the lepton mixing angles is negligible 34J. 
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TV X J 1 

Model 


10 


120 


126 


2 

Xk 


A 


3 


— 


3 


10 6 


B 


— 


3 


3 


46 


C 


3 




1 © V © 1" 


46 


D 




3 


1 © V © 1" 


46 


E' 


1 




3 


7 x 10 4 


E 


1 © V © 1" 




3 


0.89 



Table 2: A variety of renormalizable GUT models based on 50(10) x A±. Each line 
corresponds to a distinctive model scenario. Columns 2-4 specify the transformation 
properties of the 5O(10) scalar multiplets under the flavor symmetry group v4 4 . The last 
column gives the best-fit values Xd£ when fitting charged-lepton and down-type quark 
masses. 

x for the set the 14 model parameters and -P«(x) for the resulting model predictions, one 
can define a x 2 function by 

X 2 (x) = f (^^) 2 - (16) 

i=l ^ 1 ' 

The global minimum of x 2 wm represent the best possible agreement of theoretical pre- 
dictions and experimental data. This minimization task is performed using the downhill 
simplex method. 

For investigating the variation of x 2 as a function of the value O of an observable 
O, we add the "pinning term" (-P(x) — O) 2 /(0.01O) 2 to x 2 , where -P(x) represents the 
theoretical prediction for O. Note that if O agrees with one of observables Oi occurring 
in x 2 °f Eq. (TIB"]) , the Oj term has to be removed from Eq. (1161) . The small error in the 
denominator of the "pinning term" guarantees to pin the observable O down to the value 
O. The pinning procedure performs as desired when the contribution of the pinning term 
to x 2 is negligible. 

As mentioned at the end of Section[3l we have not only investigated the scenario defined 
by Eq. ([9]) but also a variety of others which are characterized by the transformation 
properties of their scalar 50(10) multiplets in columns 2-4 of Table [2] (models A-E'). We 
have found that all these scenarios fail already to reproduce the down-quark and charged- 
lepton masses — see the value of the corresponding x\t m ^ ne ^ as ^ column of Table [20 For 
comparison we have also presented the x\i °f our successful scenario in the line labeled 
by E, which will be investigated in the rest of this paper. 

5 For case A this failure is trivial: Md and Mi are symmetric with a vanishing diagonal, therefore, 
Eq. (J21) holds, which is in contradiction to the strong hierarchy in the down-quark and charged-lepton 
masses. 
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4.1 Predictions for the case of normal neutrino mass ordering 

We search for the best-fit solution for the normal neutrino mass spectrum m\ < m 2 < m^. 
In this case we find an excellent fit with the following properties: 



(17) 



X 2 = 0.96, 

m 1 = 2.795 x 1(T 2 eV, m 2 = 2.933 x 1(T 2 eV, m 3 = 5.728 x 1(T 2 eV. 

The corresponding values of the matrix elements of H, F, M. v are given by 

14.8905 
H I 14.9798 

1189.49 t 

4.45699 e i0 - 990889 - 83.3159 e^ - 940907 *" 

4.45699 e ia " 0889 " 86.5511 e ia9465577r 

83.3159 e -i °-940907t 86.5511 e ia946557 " 

2.83435 e ia26795Chr 2.93292 e ia6495677r 

M„ = ( 2.83435 e ia2679507r 2.82486 e i0 ^ \ x 10" 

2.93292 e i0M9567n 2.82486 e <aftr 

where the numerical values in H and F are in units of MeV, while the entries in A4 U are 
in units of eV. 

The non-zero value of x 2 stems from the deviation of the bottom-quark mass m& from 
its central value by +0.98er. The remaining observables of Table [1] are fitted perfectly. 
Thus the model succeeds in correcting the too large value for the mixing angle #13 of 
the Zee-Wolfenstein model, despite the close relationship between and Mi given by 
Eq. (fT3|) which, on the other hand, leads to the desired unification of m& and m T , as will 
be discussed in Section 14.31 

As explained in Section El the three light neutrino masses cannot be independent 
of each other. The sum of the eigenvalues of M. v must be zero, which translates into 
mi +1712 — 1^3 = 0. This can easily be verified for the neutrino masses of the best-fit (jPTI) . 
The sum of the neutrino masses is E = m, = 2m 3 = 0.11 eV, which lies safely below 
the cosmo logical bound E < 1 e V [35] . 

The neutrino mass spectrum has to fulfill the approximate relation ([6]). Inserting the 
central value for Am 2 tm from Table [1] into Eq. ([6]) gives m\ ~ ui2 — 2.89 x 10 -2 eV and 
7713 ~ 5.77 x 10~ 2 eV, which is in good agreement with the above best-fit results. 



The quantity R = rri\/ y Auiq measures how hierarchical a neutrino mass spectrum 
is. x 2 as a function of R is depicted in the right panel of Figure [U We read off that 
R ~ 3.1 is preferred and for the values 2.4 < R < 3.7 one obtains fits with x 2 ~ 15. Thus 
the mass spectrum is neither hierarchical nor quasi-degenerate@ but is located between 
these extrema. The narrow range of allowed values for R reflects the clear-cut prediction 
of the Zee-Wolfenstein mass matrix for the neutrino mass spectrum. 

Figure [1] (left panel) shows the constraints on the atmospheric mixing angle 823- One 
can see that values of sin 2 623 smaller than 0.38 (—2a) are strongly disfavored and thus 



^Typically, quasi-degenerate neutrino spectra would correspond to R > J Am^/Anig ~ 5.6 . 
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a strict lower bound for # 2 3 is established. However, very good fits are also possible for 
values of sin 2 #23 significantly larger than the best-fit value of 0.5. 

Concerning the solar mixing angle, Figure [1] (middle panel) shows that the whole 
physically allowed range for sin 2 9 12 gives excellent fits and therefore no prediction can be 
obtained. 

Regarding the mixing angle #13, the best-fit solution gives a value of sin 2 #13 = 2 x 10~ 4 . 
However, also significantly smaller (down to 10 -6 ) and larger values (up to 0.1) for sin 2 # 13 
are equally allowed. Thus the severe problem of the original Zee-Wolfenstein mass matrix 
(sin 2 #13 ~ 1/3) can be resolved completely by contributions from the charged-lepton 
sector. 

The best-fit gives 5pmns = 31° for the leptonic CP phase. However, varying 5pmns 
shows that the whole [0°, 360°] range allows for very good fits and therefore no prediction 
can be made. 

The effective Majorana mass of neutrinoless double-/? decay \(mpp)\ for the normal 
spectrum is given by 



where (3\ and fa are Majorana phases. Here and in the following we use the abbreviations 
C12 = cos #12, S12 = sin6 l 1 2, etc. After inserting into Eq. (fl9|) m 1 from the best-fit (TTTj) . 
employing for the other parameters the corresponding central values from Table [TJ and 
varying the two phases (5\ and (3 2 freely between 0° and 360°, we obtain the bounds 



On the other hand, the phases (3\ and /?2 are actually functions of the parameters of 
our scenario and are determined by the fit. Using the best-fit parameters ([TBI for the 
calculation of the effective Majorana mass, we obtain \(mpp)\ = 28.4 meV, which is 
identical to the upper bound in fT20l . Figure [2] presents the change of x 2 when \(mpp)\ 
is varied. We can read off that the range for the effective mass is much more restricted 
than ({IDD would suggest. Obviously, the increase of x 2 for larger values of | (mpp) | is caused 
by exceeding the upper bound of ( l20l) . The strong increase of x 2 f° r smaller values of 
\{ m 0p)\) however, is a clear-cut model prediction. For instance, allowing for only moderate 
good fits with x 2 ~ 5 results in the severely restricted range 25meV < \ (nipp) | ^ 31 meV, 
which could be tested by future neutrinoless double-/? decay experiments sensitive to 
\{m pp )\ > lOmeV. 

4.2 Predictions for the case of inverted neutrino mass ordering 

The best-fit solution for the inverted neutrino mass spectrum m 3 < m\ < m 2 turns out 
to be excellent as well. It is characterized by the following properties: 

X 2 = 0-92, , s 

(21) 

mi = 4.921 x 10~ 2 eV, m 2 = 5.000 x 10" 2 eV, m 3 = 7.963 x 10" 4 eV, 





(19) 



10.2 meV < | (mpp) | < 28.4 meV. 



(20) 
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Figure 1: 
m m 



as a function of sin 2 #23 (left panel), sin 2 6*12 (middle panel) and of R = 
Amg (right panel), where m min = mi for the normal and m 3 for the inverted 



r 



neutrino mass spectrum. The solid lines correspond to normal neutrino mass ordering, 
while the dashed lines refer to the inverted neutrino mass spectrum. 



with the matrices 
H 



My 







23.3969 


6.34572 e i7T 


116.271 

4.32801 


0.0934209 e i 



39.2455 
116.271 


2.42291 e i7r 
0.0934209 e i7r 




'.22) 



x 10" 



where the numerical values in H and F are in units of MeV, while the entries in Ai u are 
in units of eV. 

The x 2 analysis reveals that the removal of the non-trivial complex phases from F and 
M. v does not affect the goodness of the fit. Thus we specified here the fitting parameters 
for the CP conserving caseQ However, the subsequent numerical analysis is performed 
with the inclusion of the five phase parameters (CP non-conservation). As in the case of 
normal neutrino mass ordering, the main contribution to \ 2 is caused by the bottom-quark 
mass mi,, being too large by 0.95a. All the other observables are fitted very accurately. 



7 For the normal neutrino spectrum, however, the CP conserving case results in a worse, but still 
very good fit with x 2 = 1.94. Here, the main contributions to \ 2 stem from m& (+1.23ct) and sin 2 6*23 
(-0.64ct) . 
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Figure 2: x 2 as a function of the effective Majorana mass \ (mpp)\ probed in neutrinoless 
double-/? decay experiments. The solid line corresponds to normal neutrino mass ordering, 
while the dashed line refers to the inverted neutrino mass spectrum. 

Thus the GUT model allows for a considerably reduction of the maximal solar mixing 
angle #12, which spoiled the Zee-Wolfenstein model. 

M. v being of the Zee-Wolfenstein form implies m\ — m2 + 7713 = for the three light 
neutrino masses. Taking the neutrino masses from the best-fit ( 12T1) . we get £ = Xlj m i = 
2m 2 = 0.10 eV, which is safely below the cosmological limit |35j . Inserting the central 
values for the mass-squared differences from Table [1] into Eqs. (jl]) gives mi ~ m 2 — 
5 x 10~ 2 eV and m 3 ~ 7.9 x 10~ 4 eV, which is in good agreement with the numerically 
obtained best-fit values (I2TT) . 

X 2 as a function of R = m^j \J 'Am|> is shown in Figure [1] (right panel). We can 
read off that R ~ 0.09 is preferred and for the values 0.077 ^ R < 0.11 one gets fits 
with x 2 ~ 15. As in the case of normal neutrino mass ordering, the range for R is very 
restricted. Hierarchy is strongly preferred, however, too small values for m 3 are strictly 
forbidden. 

Figure [1] (middle panel) depicts the constraints on the solar mixing angle 812. We can 
read off that values for sin 2 #12 smaller than 0.3 become increasingly disfavored. However, 
very good fits can also be found for values of sin 2 6*12 larger than the best-fit value, and 
maximal solar mixing also represents a very good fit. 

Regarding the atmospheric mixing angle 623, Figure [1] (left panel) reveals that the 
whole physically allowed range for sin 2 623 gives very good fits and therefore no prediction 
can be obtained. This property is seemingly a legacy of the original Zee-Wolfenstein 
model, where the atmospheric mixing angle for inverted neutrino mass ordering is uncon- 
strained [2D] . 

As for the mixing angle #13, we find sin 2 6 13 = 2.5 x 10 -3 for the best fit. However, 
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pinning sin 2 # 13 in x 2 shows that also smaller (down to 10~ 6 ) and larger values (up to 0.1) 
are possible. For instance, enforcing sin 2 6 Vi = 10 -6 still allows x 2 = 3.1. 

Concerning the leptonic CP phase <5pmns5 the specified best-fit, which employs only 
trivial complex phases, gives <5pmns — 180°. However, varying <5pmns in the range 
(90°, 270°) allows for fits of equally good quality. Only the neighborhood of <5pmns — 0° 
seems to be slightly disfavored by x 2 — 3.2. 

The effective Majorana mass for the neutrinoless double-/? decay is given by 



With m 3 from the best-fit and taking for the other parameters in (123]) the corresponding 
central values in Table [U free variation of the two complex phases results in the following 
bounds on the effective Majorana neutrino mass: 



On the other hand, employing the best-fit parameters (T2"2"|) . we obtain \(rripp)\ = 18.4meV, 
which is located close to the lower bound of ([21]). 

Figure [2] shows the change of x 2 under variations of \(rnpp)\. We can see that the 
range for the effective mass is less restricted than in the case of normal neutrino spectrum. 
Clearly, the strong increase of x 2 f° r smaller values of \{mpp)\ comes from falling below 
the lower bound of (I24p . The rise of x 2 is l ess dramatic when moving to larger values 
of \(mpp)\. However, there is a clear bias towards values of \{rripp)\ in the lower half of 
the range spanned by fl24l) . We can also read off from Figure [2] that allowing moderately 
good fits with x 2 ~ 5 gives 13 meV < \{mpp)\ < 35 meV. Moreover, we can see that the 
\( m f3p)\ regions where x 2 ~ 2 are overlapping for both neutrino mass orderings and thus 
one cannot discriminate with \(rnpp)\ between normal and inverted mass spectrum in the 
overlap region. However, ((m^)! ^ 20 meV (which is preferred) or \{mpp) \ > 33 meV is 
only possible for an inverted hierarchy in our scenario. 

4.3 b — r unification 

As has already been noticed in Sections 14. II and 14.21 both best-fit values of mj, are located 
near the upper la bound of its input value from Table [TJ The ratio mb/m T , using the 
mean values from Table [IJ is 0.82. Employing the best-fit values for rrib and m T , this ratio 
is higher, namely mb/m T = 0.93 for both normal and inverted neutrino mass ordering. 

Figure [3] depicts x 2 as a function of rrib for both neutrino mass orderings. This figure 
clearly reflects the feature mentioned above since for values of m;, below 1190 MeV, x 2 
increases dramatically and lower values of rrib become strictly ruled out. 

There also exists an upper bound on rrib in Figure [3] at about 1250 MeV, which is 
located below the central input value of m T at 1292 MeV. In contrast to the normal 
neutrino mass spectrum, however, the inverted spectrum seems to prefer values for rrib 
near its lower bound, as can be read off from Figure [3 





(23) 



18.5 meV < \{m pp )\ < 49.5 meV. 



(24) 
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Figure 3: x 2 as a function of the bottom-quark mass mj. The solid line corresponds to 
normal neutrino mass ordering, while the dashed line refers to the inverted neutrino mass 
spectrum. The shaded region indicates the la interval for from Table [TJ 

In summary, our scenario imposes rather rigid constraints on mb and favors b — t 
unification. This feature is apparently caused by the 5*0(10) relation ([TBI between the 
mass matrices of charged-leptons and down-quarks, which differ only by a factor of —3 in 
the off-diagonal matrix elements. 



5 Conclusions 

In this paper we have presented an attempt to combine an 50(10) GUT with the family 
symmetry A 4 . We have considered renormalizable Yukawa interactions, therefore, the 
choice of scalar 5O(10)-plets for fermion mass matrices is confined to 10, 120 and 126. 
The three fermion families are accommodated in an A 4 triplet. For fitting purposes we 
use the fermion masses at the GUT scale evolved by the renormalization group equations 
of the MSSM. As a further important prerequisite we assume that the VEVs occurring in 
the scalar 5O(10) multiplets can be freely chosen for the purpose of fitting fermion masses 
and mixings. Our investigation consists of two steps — for the details see Section [3j 

In the first step we have considered only the down-quark and charged-lepton mass 
matrices. We have assigned all possible v4 4 representations to the 10, 120 and 126 and 
checked, if the down-quark and charged-lepton masses can correctly be reproduced. In 
this way we have identified a unique successful scenario given by the 10 in 1 © 1' © 1" 
of A± and the 126 in the 3 of A±. This is a non-trivial result, because we use only six 
masses for probing mass matrices constructed with more than six parameters. The mass 
matrices (1131) of the successful scenario ([9]) reflect the 5O(10) x A A structure: The 10-plets 
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contribute a general diagonal and the 126-plets a general off-diagonal matrix to the mass 
matrices. 

In the second step we have assumed type II dominance in the seesaw mechanism 
generating light neutrino masses. Since only the 126-plets contribute to the neutrino 
mass matrix, we obtain the Zee-Wolfenstein mass matrix. 

The scenario gives an excellent fit to all known data on fermion masses and mixings 
and shows, therefore, the compatibility of the family group A 4 with SO (10) GUTs. In 
summary, we have found the following features: 

• All mass matrices are symmetric. 

• In the charged-fermion sector, as an effect of SO (10) x A4, the building blocks of 
the mass matrices are general diagonal and off-diagonal matrices, generated by the 
VEVs of scalar 10-plets and 126-plets, respectively. 

• M. v is given by the Zee-Wolfenstein matrix with its definite predictions for the 
neutrino masses derived from Eq. (J2]). 

• The scenario can equally well accommodate normal and inverted neutrino mass 
spectra. 

• The lepton mixing angles of the Zee-Wolfenstein mass matrix which are in dis- 
agreement with the data are corrected by contributions to the PMNS matrix from 
M t . 

• There are definite predictions for \(mpp) \ for both spectra. 

• Our scenario gives b-r Yukawa unification. 

On the negative side we note that our scenario does not feature tri-bimaximal lepton mix- 
ing, however, some minor constraints on the lepton mixing angles exist — see Sections 14.11 
and 14.21 Moreover, up-quark masses and the CKM matrix are completely free and can 
thus be adapted to the data without imposing any restrictions on the parameters of the 
mass matrices Md, and Ai u . 

Acknowledgments: We thank L. Lavoura for valuable suggestions and reading the 
manuscript. 
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